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Abstract 

^sij \ A technique is presented to evaluate the performance of a wireless link with a multi-antenna linear Minimum- 

Mean-Square Error (MMSE) receiver in the presence of interferers distributed according to non-homogeneous Poisson 
\ processes or Poisson cluster processes on the plane. The Cumulative Distribution Function (CDF) of the Signal- 

to-Interference-plus-Noise Ratio (SINR) of a representative link is derived for both types of networks assuming 
independent Rayleigh fading between antennas. Several representative spatial node distributions are considered, for 
• which the derived CDFs are verified by numerical simulations. In addition, for non-homogeneous Poisson networks, 

it is shown that the Signal-to-Interference Ratio (SIR) converges to a deterministic non-zero value if the number 
of antennas at the representative receiver increases linearly with the nominal interferer density. This indicates that 
to the extent that the system assumptions hold, it is possible to scale such networks by increasing the number of 
receiver antennas linearly with user density. The results presented here are useful in characterizing the performance 
of multiantenna wireless networks with non-homogenous spatial node distributions and networks with clusters of 
■ users which often arise in practice, but for which few results are available. 
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Multiantenna systems have the potential to increase data rates in wireless networks through spatial multiplexing, 
beamforming and interference mitigation. The ability of antenna arrays to suppress interference is highly dependent 
psj 1 on relative signal strengths which in turn are dependent on the spatial separations between nodes. Hence, the 
^ [ spatial distribution of nodes is an important factor that influences data rates, and stochastic geometry has become 
J> ■ an important tool to analyze the performance of wireless networks as has been observed in [1] and references 
^ " therein. 

5_( \ The majority of the results in the literature that explicitly model the spatial distribution of multiantenna nodes 



have focused on homogenous Poisson spatial node distributions, i.e. systems where node positions are independent 
of one another and are distributed uniformly randomly on a plane. For instance, 13, fSj, S, S, iQ, ||7i have 
analyzed multiantenna systems in homogenous Poisson networks under a variety of assumptions using several 
different approaches. 

While simpler and (somewhat) analytically tractable, homogeneous spatial node distributions may not apply 
in many scenarios. For instance, in networks with hotspots or clusters of users, which are often encountered in 
practice, the homogenous Poisson model fails. There are also other scenarios where a non-homogenous spatial node 
distribution is appropriate such as in networks with physical limitations on user locations, e.g. roadways. 

A number of works such as |j8|, 111, fTOl . and ITTII have considered interference modeling in non-homogeneous 
Poisson and clustered networks, but for only single antenna systems. Relatively fewer works have analyzed multi- 
antenna systems in non-homogeneous or clustered networks. |12| considered interference-alignment in clustered 
wireless networks where a partial interference-alignment protocol is used to reduce the system to a form similar 
to a single-antenna system. While interference-alignment can provide enormous data rates, it requires significant 
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overhead for the exchange of transmit (Tx) Channel-State Information (CSI). In comparison, the system we analyze 
is more attractive for implementation as it does not require Tx CSI and uses a linear receiver which only requires CSI 
of the target transmitter and the spatial covariance matrix of the interference plus noise. The spatial interference 
plus noise covariance matrix can be approximated by a sample covariance matrix which can be obtained by 
appropriately averaging aggregate received signals. |[T3ll approximates networks of multiple-input-multiple-output 
(MIMO) links with Carrier-Sensing-Multiple-Access (CSMA) using a Poisson approximation for the spatial node 
distribution. However, the multiple antennas are not used for interference mitigation compared to our work, which 
does. 1 14] and [15] respectively consider non-homogenous Poisson networks and networks with certain forms of 
spatial correlation of node positions, in particular hard-core networks. Although 1141 overlaps considerably with 
this work in the system model, 1 14] uses an asymptotic analysis which is applicable only with moderately large 
numbers of antennas as compared to our work which applies to all numbers of antennas. Furthermore, they do 
not provide the distribution of the SINR, instead focusing on the convergence of appropriately normalized versions 
of the SINR as the numbers of antennas per receiver gets large. In the equivalent asymptotic regimes, our results 
agree with those findings. 

In this paper, we develop a framework to characterize the SINR of a representative link between a multiple-antenna 
receiver and a single antenna transmitter in the presence of single-antenna interferers that are distributed according 
to a non-homogeneous Poisson process on a plane, whereby node locations of the interferers are independent of 
one another but the spatial distribution of these nodes is non-uniform (e.g. the node distribution in Figure [T]). We 
assume independent Rayleigh fading channels among all antennas and use the linear MMSE receiver which is the 
optimal receiver (to maximize SINR) amongst all linear receivers. These results are then extended to networks 
where the non-homogeneity of spatial node distributions of the interferers is random, and in particular, to Poisson 
cluster processes. Note that this latter model is no longer Poisson, as node locations of the interferers could be 
correlated. For instance in clustered networks with random clusters, the presence of a node at a particular position 
increases the likelihood that there are other points nearby. In our analysis, we shall focus on two types of Poisson 
cluster networks, namely the Matern and the Thomas cluster processes which have been analyzed for single antenna 
systems in |9], although our methods could be extended to other models as well. 

For non-homogeneous Poisson networks, the node distribution of the interferers is modeled by a deterministic 
intensity function A(r, 9) which captures the likelihood of interferers occurring in an infinitesimal region around 
a point {r,9). The intensity function is used to develop a framework to compute the CDF of the SINR of a 
representative link in a network with interferer distribution governed by A{r,9). With this framework, we find a 
closed-form expression for the CDF of the SINR of a representative link in the center of a cluster with a power-law 
distribution of interfering node intensities (analyzed asymptotically in [14]). A particular example where A(r, 9) = 
0.1r~^ is shown in Figure[T] We use a high density here to emphasize the clustering effect at the origin. Furthermore, 
we also provide expressions involving generalized functions for networks with circular-symmetric Gaussian intensity 
functions, piecewise power law intensity functions, and intensity functions expressed by polynomials within a 
bounded region. This latter result can be used to approximate with arbitrary accuracy, a large class of continuous 
intensity functions, by using polynomial approximations. 

The non-homogenous Poisson results are extended to systems with interferers distributed according to a Poisson 
cluster process by first conditioning on one realization of the cluster centers, which results in a non-homogeneous 
Poisson process of interferers. The framework developed for non-homogeneous Poisson processes is applied, and 
the conditioning is then removed to obtain the unconditioned CDF of the SIR on a representative link. Note that 
we assume the noise is negligible for the Poisson cluster model to simplify analysis. 

Figure [2] illustrates a Matern cluster network where clustered interferers (points) surround a representative receiver 
(triangle). The crosses represent cluster centers which follow a homogeneous Poisson process, and interferers are 
uniformly distributed in a disc of a fixed radius in each cluster. Please disregard the labels for Case 1 and 2 in the 
figure for now as they will only be used later in the paper. 

Additionally, using the CDFs we derived for non-homogenous Poisson networks we show that as the number 
of receiver antennas L is scaled linearly with A(r, 6) in non-homogenous Poisson networks, the SINR converges 
in probability to a positive deterministic value which is a function of A{r,6). This result indicates that to the 
extent that our assumptions (in particular independent Rayleigh fading) hold, increasing the number of antennas 
per receiver can help scale such networks by supporting increasing user densities. This finding generalizes similar 
findings derived previously for homogenous Poisson networks in 131, IH and IS. 



Fig. 1. Randomly distributed interferers (points) with intensity function A(r, 9) = O.lr ^ , with the representative receiver (triangle) at the 
center of the network. 



Fig. 2. Homogeneously distributed interferers (points) restricted to a disc at each of the homogeneously distributed cluster centers (crosses) 
with the representative receiver (triangle) at the center of the network. Case 1 represents a cluster that does not include the origin, whereas 
Case 2 represents a cluster that includes the origin. 

The results in this work are derived combining techniques developed for homogeneous networks in ||6l and non- 
homogeneous Poisson network models from |14|. The extensions to Poisson cluster networks is made by applying 
techniques from stochastic geometry including an apparently novel form of Campbell's theorem (e.g. see |16|). 

Section |ll] introduces the system model with Section III-AI describing the general system model shared by both 
the non-homogenous Poisson and Poisson cluster networks. Sections IlI-BI and IlI-CI contain elements of the system 
model specific to non-homogenous Poisson networks, and Poisson cluster networks respectively. Sections Hill and [TVl 
provide the main results for non-homogenous Poisson networks, and Poisson cluster networks respectively. Section 
|V] presents applications to representative intensity functions; Section [Vlll illustrates numerical results for different 
specific models. Finally, Section IVIIII summarizes and concludes this paper. Throughout the paper, uppercase bold 
characters represent matrices and lowercase bold characters represent vectors. The indicator function equals 
1 if ^ is true, and otherwise. 



A. General Channel Model 

In this subsection, we present elements of the system model that are shared between the non-homogenous 
Poisson and Poisson cluster models. A representative receiver located at the origin is communicating with a 




II. System Model 



representative transmitter at a fixed distance rj-. Note that since the spatial distribution of interferers in the network 
is non-homogenous, the representative receiver does not correspond directly to the notion of the "typical" receiver 
commonly encountered in the literature (e.g. see [1] and references therein) because in non-homogenous networks, 
statistical properties of the system at any point on the plane (e.g. the origin) could differ from the properties at 
other points. For the purposes of this work, the representative receiver should be interpreted simply as the receiver 
at the origin, and the representative transmitter is the transmitter to which it is linked. 

The representative link is surrounded by interfering nodes distributed on a plane, transmitting simultaneously in 
the same frequency band. We assume that these interferers, with equal transmit power, are communicating with 
other receivers whose locations do not affect our results. Loss in signal power due to propagation is modeled by 
the inverse power-law model, so that the average power (over fading realizations) p from a node transmitting with 
unit power, received at a distance r is p = r^", with the path-loss exponent a > 2. The receiver has L antennas, 
and the representative transmitter and each interferer has a single antenna. We use the subscript T to denote the 
representative transmitter and 1,2,- •• ,n to label the interferers. Using this notation, rj represents the distance 
between the i-th interferer and the representative receiver at the origin. The transmitted scalar symbol from the 
representative transmitter is xt and the symbol from the i-th interferer is Xj. The L x 1 received signal vector can 
be described as follows: 

n 

-a/2 , ~a/2 , 

y = rj. ' gTXT + y gia^i+w, (1) 

i=l 

where Tj^'^^'^^t (or r- °''^gi) represents the channel coefficients between the representative transmitter (or the i-th 
interferer) and the receiver, g^, gj G C^^^ consist of independent and identically distributed (i.i.d.), zero-mean, 
unit-variance complex Gaussian entries, vi^ is a noise vector with zero-mean, i.i.d. complex Gaussian entries with 
variance per complex dimension. 

The representative receiver estimates xx from y using a linear MMSE estimator that is known to maximize the 
SINR. The SINR is given by: 

SINR = r^^gj. (cPGt + a^Ii) g^, (2) 

where I^, is an L x L identity matrix, and P = diag [rj"", r^", • • • r~"] is a n x n diagonal matrix with entries 
corresponding to the path losses from the interferers to the representative receiver. The i-th column of G G C^^" 
is gj. To simplify notation, we define the distance-normaUzed SINR as 7 = SINR • r^. 

We consider this same general channel model for two different spatially distributed interference models: non- 
homogeneous Poisson models and Poisson cluster models. 



B. 1^ on-homogeneous Poisson Networks 

For the non-homogeneous Poisson models, the interferers are distributed according to a non-homogenous Poisson 
point process (PPP) on the plane with intensity function A(r, 6), which is expressed in polar coordinates. In order to 
construct the non-homogenous PPP of interferers, we start with a circular network of radius R and place interferers 
at i.i.d. random locations with the probability density function (PDF) frfi{f,0) which is related to the intensity 
function as follows: 

The number of interferers n in the circle with radius i? is a Poisson random variable with mean /i, expressed as: 

f-R 1-271 

fi= / rA{r,e)dOdr. (4) 

In the derivation of the main results we take i? — )• 00 to model the interferers distributed according to a non- 
homogeneous PPP with the intensity function A{r,9). 



C. Poisson Cluster Networks 

To model networks with random clustering effects, we assume interferers surrounding the representative receiver at 
the origin are distributed according to a Cox process (also known as a doubly stochastic process), which generalizes 
the non-homogenous Poisson process in that the intensity function is random (e.g. see [17]). A particular realization 
of the intensity function of the Cox process results in a non-homogenous Poisson point process. In particular, we 
consider Neyman-Scott processes, which are Cox processes with clusters whose centers are distributed according 
to a stationary PPP (see e.g. fTTl). Note that for simplicity, the representative transmitter is not assumed to belong 
to any cluster in our model and is fixed at a distance tt from the representative receiver. Additionally, we shall 
assume that the noise power is negligible compared to the interference in our model, which thus applies to networks 
with a high density of interferers. This assumption is made to simplify the analysis. 

To define the Neyman-Scott process within our context, we start with the intensity function associated with a 
single cluster whose parent point (or cluster center) is Xi, defined as A(r, ^; Xj). The parent point is surrounded by 
i.i.d. daughter points, the number of which is a Poisson random variable with mean /z^. The relationship between 
the PDF of the points associated with a cluster at Xi and the intensity function of each cluster depends on the 
clustering model and is given in Section |Vl] for the Thomas and Matern cluster processes considered in this paper. 

Alternatively, the daughter points can be constructed by first generating them around the origin and then displacing 
the resulting points by the associated parent point. In other words, to generate the Ui daughter points associated 
with cluster i, which we denote in Cartesian coordinates by Yii,Yi2 ■ ■ -Yin^, we first generate a Poisson random 
variable n-i with mean /i^. Then, we generate rij points l^i, 1^2 • • • Yirn i-i-d. with PDF 

/r,e(r,0) = — A(r,^;0) (5) 
fJ-d 

which are essentially daughter points with a parent point at the origin. These points are then displaced by Xi to 
produce the true daughter points, i.e. Yij = Yij + Xi. 

The set of parent points, denoted by 11 = {Xi\ i = 1,2, ...}, is generated from a homogeneous PPP with intensity 
Pp on the plane. All the daughter points are considered interferers in the network, but the parent points are not. 

Thus, the intensity function the interferers in the network, conditioned on a particular realization of the parent 
point process 11, is: 

A{r,e;U) = Hr,0;Xi). (6) 

When the conditioning on the realization of the parent point process is removed we obtain the Neyman-Scott cluster 
process. 

To illustrate, consider Figure |2] which illustrates a realization of the Matem cluster process (a type of Neyman- 
Scott process) where each cluster is simply a disk of radius R^. The corresponding per-cluster intensity function 
is 

A(r, 6; Xi) = Pd'i-{{r,e)eB{x,,Ra)} > (7) 

where B{Xi, R^) denotes a radius R^ disk centered at Xi and pd = pd/ (^^1) is the density of the daughter points 
in a single cluster. Conditioned on a particular realization of the parent points, the Matern cluster process is thus 
a non-homogenous PPP with intensity function: 

Air,e;U)= pdl{(^ 

III. Main Results for Non-homogeneous Poisson Networks 

In this section, we consider non-homogenous Poisson networks described by deterministic intensity functions 
A(r, 9). We first derive the outage probability of the SINR under assumptions in lll-Al and lll-Bl in our system model. 
Then, with the same system model, we show that by scaling the number of receiver antennas linearly with the 
intensity function, the SINR on the representative link converges in probability to a deterministic value determined 
by the "shape" of the underlying intensity function. 



A. Outage Probability 

One of the key performance measures of communication systems is the outage probability, which is typically 
defined as the probability that the SINR is below a threshold r. For a fixed r^, this probability is PrjSINR < r} = 
Fj{Tr^), where ^7(7) is the CDF of the distance-normalized SINR 7 given by the following theorem. 

Theorem 1: The CDF of 7 in a network with non-homogeneous Poisson distributed interferers is 

L-l 



^ r(L,^(7)+a27) 



r(L) 

where 



(9) 



V^(7) = y^ Air,9)r ^^^_l^ d9dr, (10) 

and r(.) and r(., .) are the gamma function and the upper incomplete gamma function. In addition, the corresponding 
PDF of 7 is: 

(^(7) + a^7)^~^ exp(-V^(7) - a^j)ia^ + -^^7)) 
Ml) = ^3T)^ , 

where ^'(7) is the derivative of ^"(7) with respect to 7. 
Proof: Given in Appendix [A] 

With this result, we can quantify the effect of having more antennas in non-homogeneous networks. For instance, 
increasing the number of antennas at the receiver from L to L + 1 reduces the outage probability by (^'(7) + 
(T^7)^exp(--(/'(7) - CT^7)/L!. 



B. Scaling Non-homogeneous Networks by Increasing the Number of Antennas 

One of the questions that Theorem 1 allows us to answer is whether one can maintain a non-zero SINR if the 
number of antennas at the representative receiver is increased linearly with the nominal density of interferers in 
the network. In the context of homogeneous networks |3|, |5| and [6| found that this is indeed the case. Here, we 
shall show that a similar result holds even when the spatial interferer distribution is non-homogenous. 

Assuming that noise is negligible, we show that the SIR on the representative link converges to a deterministic 
value if the number of antennas at the receiver increases linearly with a nominal interferer density. This is under the 
assumption that the channel model, independent Rayleigh fading in particular, holds, and that accurate measurements 
of CSI are available at the receiver. A key result that we use is the following lemma which may already be known 
but we were unable to find it in the literature. 

Lemma 1: Let the upper regularized gamma function be denoted by Q{L,x) = ^f^y, where T{L,x) is the 
upper incomplete gamma function. Let L be a positive integer and g > 0, then 

lim Q(L,gL) = (°' ^ " J (11) 
L^oo I 1, if g < 1. 

Proof: Given in Appendix |B] Note that the proof here is a coiTccted version of the proof of Lemma 1 in a 
conference version of this paper, fTSl. 

Suppose that the intensity function K{r,9) = f3Ac{r,6), where Ac{r,6) is a nominal intensity function which 
describes the "shape" of the true intensity function, and /3 is the nominal interferer density which scales the nominal 
intensity function. We also define: 

Ml) = ^c{r,e)r ^^^_l^ dedr. (12) 

Theorem 2: Let (5 = iL with a constant scaling coefficient £ > 0. As L — ;> 00, the SIR converges in probability 
to i;-^ (i) r--. 

Proof: Given in Appendix O 



Therefore, if we increase the number of antennas Unearly with the nominal interferer density in a non-homogenous 
Poisson network, the SIR will approach a constant non-zero value. This fact implies that such networks can be 
scaled by linearly increasing the number of antennas per receiver with user density without degrading the SIR 
to zero, provided that the assumptions of the system are satisfied. Note that as the number of antennas gets very 
large, the independent Rayleigh fading and accurate channel state information assumptions will become increasingly 
difficult to satisfy as they would require increasing antenna separations and increasing channel estimation times. 



IV. Outage Probability for Neyman-Scott Cluster Networks 

In this section, we derive the outage probability of the SIR in a network with interferers distributed as a Neyman- 
Scott cluster process under the assumptions in Sections III-AI and III-CI Even though these networks are described 
by nondeterministic intensity functions due to the random locations of parent points, the Neyman-Scott process 
conditioned on one realization of the process of parent points becomes a non-homogenous PPP Thus, the approach 
we take is to first condition on a realization of the parent point process, and then apply Theorem [T] to the resulting 
non-homogenous PPP. Finally, we remove the condition by taking the expectation over the parent point process 
using the following Lemma that may already be known, but we were not able to find in the literature. 

Lemma 2: Let 11 be a Poisson process on with density p, and let H = Xlxen ^i-^)^ where C : — (0, oo) 
is a non-negative measurable function. Then, for non-negative integers k, 

.,exp|/,(e-*.-l)4 V .[J...C^(^).-«'W]'"' 
Ur2^ ^ J ^ mi!l!™im2!2!"^^..mfc!A:!™'= 

(mi,...,mfc)eA^fc 



En 



where u and z are integration variables, and is the set of all fc-tuples of nonnegative integers (mi, m/j) 
satisfying the constraint: 

1 • mi + 2 • 7712 + 3 • m3 + ... + k ■ = k. (14) 

Note that C,^{z) refers to the j-th power of the function C,{z). 
Proof: Given in Appendix |D] 

To find the CDF of the SIR for the Neyman-Scott model, set 

aXi)= / K{r,e-Xi)r \ dOdr. (15) 

Conditioned on the process of the parent points 11, ip{^) in Theorem [T] can be expressed in terms of the intensity 
function A(r, 6';n) in ^ as follows: 

V^(7;n)= / / Mr,0;X,)r—-^dedr= J] dX,). (16) 

The interchange of the order of integration and summation follows from Theorem 11.30 in |fT9l . Then, substituting 
n) for ^"(7) in ® yields the conditional CDF of 7 conditioned on 11, i.e. F^{'y\ 11). To remove the condition, 
we take the expectation with respect to H, i.e. F^i^j) = En [-^7(7! H)], which leads to the following theorem. 

Theorem 3: Neglecting the thermal noise power, the CDF of the distance-normalized SINR 7 in a network with 
interferers distributed according to a Neyman-Scott cluster process is: 



F-^(7) = 1 - exp |27r^ ppr (e"'^^'') - l) dr| 



^ mi!l!™im2!2!™-2...m^!^!mt, 

A:=0 (mi,...,mfe)eXfc 

Proof: Given in Appendix 10 

Note that equation (fTTl ) requires the evaluation of integrals of the two forms: J^2 Pp (e"^^''^ — l) du and PpC^ {z)e~'^^^^dz. 
These integrals can be evaluated numerically using standard methods for both the Matern and Thomas cluster 
processes as given in Sections |Vl] and IVII-DI 



V. Applications to Representative Non-homogeneous Poisson Networks 

In this section, we apply the results for the non-homogenous Poisson networks given in in Section |lll] to some 
representative intensity functions. We start with polynomial intensity functions within a bounded domain which can 
be used to approximate a large class of functions. 



A. Polynomial Intensity Function in a Bounded Domain 

Let us define a polynomial intensity function in a bounded domain as follows: 

m 

^mir, 0) = Y^ akrH{o<^<R^y + por'l{R„<r}, (18) 

fc=0 

where ao, oi, 02, ■■■am are arbitrary polynomial coefficients. The intensity function follows a polynomial when 
< r < Rq, and has a power-law decay when r > Rq, for e < and po ^ 0. The power-law decay is added to 
the polynomial intensity function here for generality. Then, the corresponding CDF of 7 is: 

Fj^M = 1-1: M(2i±i:!T): ^^p(_^P(^) _ 

i=0 

where 

V'^(7) = | AP{r,e)rr^^^^dedr. (20) 

This integral is evaluated in Appendix IH which yields: 

where 2-^i('5 •; S •) is Gauss's hypergeometric function (e.g. see f20'|). By the Stone -Weierstrass Theorem fT9l, any 
continuous function can be uniformly approximated with arbitrary accuracy by a polynomial in a bounded interval. 
This fact combined with (|2T]) can be used to derive the following Theorem. 
Theorem 4: For any intensity function of the form 

A(r, 6) = /i(r)l{o<r.<Ro} + Po?''l{i?o<r} > (22) 
where h{r) is any continuous function of r, there exist coefficients ao, ai, a2^ ■■■a-m such that 

lim F^(7) = F(7). (23) 

m— >oo 

Proof: Given in Appendix [G] 

Theorem |4] allows us to approximate, with arbitrary accuracy, the CDF of the SINR for any intensity function 
that is continuous in r within a finite domain using a polynomial expression. In particular, since efficient algorithms 
exist for fitting polynomials to real data, this technique could be useful to analyze networks whose geometrical 
characteristics are not easily captured by mathematical models, but for which experimental data on positions of 
the interferers are available. Additionally, this technique can be used to approximate systems where the intensity 
function is complicated. An example of such an application is given in Section I VII- A I 



B. Piecewise Power-law Intensity Function 

Consider networks where node distribution of the interferers vary in different radial ranges according to power 
law functions. As we see in the following subsections, such a model can be used to describe networks with a dense 
cluster at the origin, both for infinite and finite-sized networks. This model could also be used to analyze networks 
with guard-zones around receivers, which will be explored in future work. 



Consider a set of non-negative numbers representing radial ranges, Rq < < 
intensity function has the following form: 



< Rm- Assume that the 



A{r,e) 



for Ro<r<Ri 
for Ri <r < R2 

for Rrn-l <r <Rr, 



(24) 



where ei > —2 if Rq = 0, and < a — 2 if R^ = 00. In the range Rk-i < r < R^, the intensity function of the 
interferers follows a power-law distribution with nominal density pk, and exponent e^. Applying Lemma [3] from 
Appendix |F] to (|24l) . we find that the CDF of 7 is given by Q with 



E 



fc=i 



2 + efc 



Rl'^\Fi 



ek + 2 ek + 2 + a RZ 



a 



efc + 2 ek + 2 + a 



a 

na 



7 



a a 7 

For the simplest scenario, we have only one piece for the intensity function as follows: 

K{r,e) = pr", 



(25) 



(26) 



where — 2 < e < 0, to maintain finite interference for any r. This intensity function can be used to model a network 
with a dense cluster of interferers, centered on the representative receiver and is useful to model networks with 
hot-spots. In this case, following the derivation in Appendix |F1 Tp{'y) in Theorem [T] can be expressed in closed form 
as: 

fOO r^TT „ 1^2, 



^"(7) 



7 



JO 

Substituting ( |27] ) into Q, we have the CDF of 7 as 

.2 



^7(^+^)/"csc(vri±^ 
a \ a 



(27) 



^^7(7) = 1 



^ il 

i=0 



CSC vr 



a 



e + 2 



a 



X exp 



271^ p f e + 2 
CSC vr 



a 



a 



^ie+2)/a _ ^2^ 



(28) 



Note that an asymptotic analysis was used to show convergence of an appropriately normalized version of the SIR 
for this power-law interferer distribution model in [14J. To confirm that the result above agrees with the conclusions 
in that work, we first neglect the noise by setting a"^ = 0. Since the SIR grows without bound as L — )• 00, we need 
to apply a normalization to the SIR as is done in 11141 to avoid degenerate results as L — )• 00. Define a normalized 
version of the SIR, ^ = L-"/(2+^)r^SIR = L-°'/^'^+''>-f. Then, 



-a/{2+e) 



7 < = Pr(7 < 



a/{2+e) 



^ r(L,g^csc(^^)e(-+^)/"L) 



r(L) 



(29) 



Given Lemma[T] if we set q = esc (vr^^) F^iO approaches a step at a deterministic value 

as L — )• 00. This implies that for large number of antennas L, SIR 
consistent with the findings in |[T4l . 



esc 



^cse(vr^] 



rj," which is 



C. Gaussian Intensity Function 

Suppose that the receiver is located in the center of a cluster whose intensity function follows a circularly- 
symmetric Gaussian function. A(r, 6) can be written in the form of the PDF of Rayleigh variable multiplied with 
some constant p expressed as: 

K{r,e)=p^e~^''^-\ (30) 

where v controls the width of the Gaussian intensity function. Then the mean number of interferers in the network 
is 

/■oo r2-n 2 

p= / p—e-'^'/'^'''dedr = pvTrV2^. (31) 
Jo Jo ^ 

As a result, the PDF of the distance from each of the interferers to the receiver can be characterized by a Maxwell 
distribution, and the PDF of (r, 9) is given by: 

2ir\ IT 

For specific values of the path-loss exponent a, the integral of V'(7) can be evaluated in terms of generalized 
functions. When a = 3, 4 and 5, we respectively have: 



fr,e{r, e) = —\l ^3 . (32) 
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(34) 
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where G represents the Meijer G-function, and 1F2 is the generalized hypergeometric function which can both be 
evaluated efficiently in most mathematical software packages. The CDFs can be found by substituting V(7) iiito 
(111). Note that for general values of a, one can use a polynomial approximation to very accurately approximate the 
CDF according to Theorem IH as illustrated in Section IVIII 



VI. Representative Neyman-Scott Cluster Processes 

In this section, we apply the techniques in Theorem [3] to two representative Neyman-Scott cluster processes, 
namely the Thomas and Matern cluster processes. These processes have been used as models for wireless networks 
with clustered interferers in works such as |9]. Recall that the Neyman-Scott cluster process uses an underlying 
Poisson point processes to model parent points of clusters, and each cluster has a random number of daughter 
points distributed according to some intensity function, centered on the parent point. 



A. Thomas Cluster Process 

For the Thomas cluster process (e.g. see ifTTl ). the daughter points are distributed in an i.i.d. fashion according 
to a two dimensional Gaussian distribution (with width parameter u) which is centered at the parent points, Xi. 
The distance between a random point from a given cluster to the representative receiver at the origin thus follows 
a Rician distribution since the points from this cluster are distributed according to a non-centralized, symmetric, 
2-dimensional Gaussian distribution. (Note that this differs from the model in Section IV-CI because there, its the 
intensity function that follows a circularly symmetric Gaussian distribution.) Hence, the PDF of the distance from 
an interferer at cluster Xi to the representative receiver at the origin can be directly written as: 

fM Xi) = ^ exp j Jo j , (36) 



where /q is the zeroth-order modified Bessel function of the first kind. From ([5]), the relationship between fr{r\ Xi) 
and A{r,9;Xi) can then be expressed as follows: 

r rA(r, 9; X,)de = fi^ Ui'^, e\ Xi)d9 = ^d/r(r| Xi). (37) 
Jo Jo 

Therefore, by the definition of C{Xi) in (fTSl ). we have 

C(Xi) = / exp —2 lo — ^ , (38) 



Substituting this function into (fTTl) yields the CDF of the distance-normalized SIR, 7. The CDF can be efficiently 
evaluated numerically as illustrated in the next section. 



B. Matern Cluster Process 

Recall the definition of the Matern cluster process in Section III-CI with the intensity function of a single cluster 
associated with parent point X^ given by ([7]). To find the corresponding PDF of the distance of a random point from 
this cluster to the representative receiver fr{r\Xi) , we need to consider two disjoint and independent cases. The 
cases correspond to whether the cluster under consideration includes the origin or not, as shown in the illustration 
in Figure [2] 

Case 1: \Xi\ > Rd, i.e. when the representative receiver is outside the disk B{Xi,Rd). In this case, r is the 
distance between a random point inside a circle of radius and a fixed point outside this circle (at a distance 
\Xi\). This PDF is given in lH as follows: 

2r /r2 + |Xi|2 - i?2\ 
fM \Xi\ > Rd) = ^ COS^^ • l{\x,\-R,<r<\X,\+R,} (39) 

Case 2: \Xi\ < Rd, i.e. when the representative receiver is inside disk B{Xi,Rd). Using elementary geometry 
and applying the techniques used to derive ( [39l) in l2TI , one can find the PDF of r in a similar form as follows: 

2r 2r /r2 + |XJ2 — i?2 \ 

fr{r\ \Xi\ <Rd) = -^^- l{r<R,-\X,\} + ^ ^ J • l{R,_|X.|<r<i?„+|X,|} • (40) 

Hence, similar to the derivation of (1381 ). substituting ( [39l ) and (l40l) into (ITSl) yields. 

C{Xi) = 



J|X,|-_Rd ' ^"'^ 1^ 2r\X,\ ) l+r—-f"'' , ^ 

Jo ^'^l+r-°7"'^ + jRi-|Xi| I 2r\Xi\ jl+r-'-y"'^ 



\XA < Rd. 



-\Xi\ ^' \ 2r\Xi\ 

The CDF is found by substituting this expression into ([TTll . The resulting CDF can be numerically evaluated 
efficiently as illustrated in the following section. 



VII. Numerical Results 
A. Monte-Carlo Simulation of Network with Gaussian Intensity 

In this section, we show numerical results for the cases analyzed in the previous sections. We ran Monte-Carlo 
simulations for different network scenarios to validate the general technique in Theorem [T] for non-homogeneous 
networks and to verify the extension of this technique in Theorem [3] for Neyman-Scott clustered networks. 

To begin, we simulated the scenario with the Gaussian intensity function described in Section IV-CI for path-loss 
exponents a = 3, 4 and 5. In addition, we used the following parameters: mean number of interferers fi = 2000 
in the network, v = 500 , the distance between the representative transmitter and receiver r-p = 20, number of 
receiver antennas L = 10, and noise variance = 10^^^. The number of interferers n was generated as a Poisson 
random variable with mean fi in each trial of the simulation, and the distances from each interferer to the receiver 
were generated according to the Maxwell distribution with parameter v. 

Figure [3] shows the simulated PDFs from three Monte-Carlo simulations with 100,000 trials each, along with 
theoretical PDFs from equations ( [33] ). ( [34l ) and ( [35l ). Additionally, we also computed a polynomial approximation 
of the PDFs from Theorem ID The approximation is based on 8 terms for the polynomial intensity function in (|2TI ). 
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Fig. 3. Comparison between the empirical and theoretical PDF of the SINR for the Gaussian intensity function for a — 3,4, and 5, with 
mean number of interferers /i = 1000, v = 500 , rr = 20, L = 10, and = 10~^*. 



Note that the theoretical PDF can be expressed for specific values of a, but the polynomial approximation does 
not require a to be known. Hence the polynomial approximation could be used in more general scenarios. The 
figure clearly validates that the simulations match the theoretical predictions, and that the polynomial approximation 
holds. 

Additionally, notice that even with a large number of trials, the simulated PDF is not smooth. This suggests that 
a purely simulation-based approach to estimate the PDF of the SINR in this case is computationally prohibitive as 
an extremely large number of trials would be necessary to produce accurate CDFs and PDFs. Hence, the CDFs 
and PDFs given in terms of generalized functions through (l33l) . (l34l ) and (l35l) . or the polynomial approximation 
from Theorem |4] are useful as they can be evaluated numerically more efficiently than running large numbers of 
simulations for each set of system parameters that we wish to characterize. 

B. Scaling interferer density by increasing number of antennas 

Section lTlI-Bl shows that if the number of antennas is scaled linearly with the intensity function, the SIR approaches 
a deterministic, positive value. With the Gaussian intensity validated above, we confirmed Theorem |2] by comparing 
the CDF of the SIR with A(r,6') = /3Ac(r,6l) = ^^e-""^ /'^''\ We consider L = 1,5,10 and 20 antennas, with 
corresponding mean number of interferers = 100, 500, 1000 and 2000 in the network. Note that these values 
correspond to a linear increase in interferer density with the number of antennas. The CDFs are illustrated in 
Figure |4l Figure |4] shows that as the number of interferers increases from 1 to 20 with a corresponding increase 
in the density of interferers, the CDF of SIR approaches a step function, i.e the SIR approaches a deterministic 
non-zero value in distribution. Moreover, the SIR converging in distribution to a deterministic value implies that it 
converges in probability. 

C. Uniform Versus Power Law Clustered Networks 

We can use the power law intensity function in ( |26l ) with different values for the exponent e to compare the SINRs 
between networks with uniform and clustered interfering node distributions with deterministic clusters centered at 
the origin. Figure |5] shows the simulated and theoretical PDFs for the intensity function A(r, 6*) = The 
additional parameters of the simulation are given in the caption. From the graph, it is clear that the theoretical 
prediction of Theorem 1 is accurate. To compare networks with different values of e, we note that if the nominal 
density p is fixed, different values of e result in radically different numbers of interferers in the vicinity of the 
representative receiver. To make a fairer comparison between networks with different degrees of clustering, we 
consider a finite, circular network of radius Rc with a fixed mean number of points, but distributed with different 
degrees of clustering. The piecewise power-law model introduced in Section IV-BI can be used here. The intensity 
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Fig. 4. Cumulative distribution function of SIR (dB) with number of antennas increasing linearly with nominal interferer density, with 
« = 500 , TT = 20, a = 3 and = 10"". 
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Fig. 5. Comparison between the empirical and theoretical probability density function of SINR with the power-law intensity function 
A{r,e) = . The parameters used are rr = 10, L = 10, a = 4, = 10"^^, and 100, 000 Monte-Carlo trials. 



function equals pr^ for < r < Rc and zero otherwise. The values of p and e are related to each other such that 
the mean number of interferers in the network is p. 

Assuming Rc = 1000 with e varying from -1 to 0, we selected values of p so that p = 3142. Figure |6] shows the 
probability that the SINR is less than or equal to 10 for different values of e and the remaining parameters given 
in the caption. Note that the outage probability increases significantly with clustering. For instance, with L = 4 
antennas, the outage probability is below 10^'^ for e = 0, while it is greater than 0.1 for e = —0.5. Additionally, 
observe that it is possible to significantly reduce the outage probability by increasing the number of antennas at the 
receiver. For instance, L = 12 antennas at the receiver in a clustered network with e ^ —0.5 has the same outage 
probability as L = 4 antennas in a homogeneous network. This indicates that increasing the number of antennas 
provides a large decrease in the CDF at the center of dense networks. 

D. Poisson Cluster Networks 

In each trial of the Monte-Carlo simulations for both Thomas and Matem cluster processes, we start with an i.i.d. 
homogeneous Poisson process with density pp for cluster centers. Then we independently generate a Poisson number 
of daughter points with mean p^ at each of the cluster centers according to the corresponding PDFs provided in 
Section |Vll Even though we could not solve the integrals in the CDFs of the SIR for either case in closed form, 




e 



Fig. 6. Outage probability vs. e for power-law intensity functions with I/ = l,2,4, 8 and 12 receiver antennas. Representative link length 
rr ~ 5, mean number of interferers in a radius 1000 circle is fixed al fi — 3142, a — A, and = 10~^^. 

we were able to apply standard quadrature numerical integration methods to evaluate, with sufficient accuracy, our 
theoretical prediction to compare with the empirical CDFs of the SIR. 

For both Thomas and Matern cluster processes shown in Figure |7] and Figure [8] respectively, we assume the 
distance between the representative transmitter and receiver is rx = 10, path-loss exponent a = 4, and zero noise 
power cr^ = to focus on interference-limited clustered networks. In the Thomas cluster process show in Figure 
|71 the homogeneous density of the cluster centers is pp = 3 x 10^^, the mean number of points at each cluster 
= 100, and the width parameter for the Gaussian density, u = 500. In the Matern cluster process shown in 
Figure m the homogeneous density of the cluster centers is pp = 1.5 x 10~^, the mean number of points at each 
cluster fiii = 200, and the radius of each cluster = 100. Then, we compare the empirical CDF and theoretical 
CDF for different number of antennas L = 1,2,3 and 4 for each case. In both cases, the simulated and theoretical 
CDFs of the SIR are virtually indistinguishable, which confirms both the accuracy of the analysis and the numerical 
integration which was performed using a simple quadrature integration. Notice that the variance of the SIR for the 
cluster processes are large due to the high degree of irregularity in the spatial interferer distribution. 




SIR(dB) 

Fig. 7. Comparison between the empirical and theoretical CDF of SIR when interferers are generated from a Thomas cluster process. 
Varying the number of antennas L, the empirical CDF results from 100, 000 Monte-Carlo trials. Other parameters used are rr = 10, a = 4, 
o-^ = 0, Pp = 3.18 X 10"^, Hd = 100 and ly = 500. 

To get a better understanding of the effects of random clustering, we compared systems with the same effective 
density of interferers but with different densities of parent points (i.e. different density of clusters). With same values 
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Fig. 8. Comparison between tlie empirical and theoretical CDF of SIR when interferers are generated from a Matern cluster process. 
Varying the number of antennas L, the empirical CDF results from 100, 000 Monte-Carlo trials. Other parameters used are rr = 10, a = 4, 
= 0, Pp = 3.18 X 10"^ Hd = 100 and Rd = 100. 
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Fig. 9. CDF of the SIR for Matern cluster processes with varying parent point density, pp and mean number of daughter points per cluster 
fid, but the same total density of 10^^. Other parameters used are rr = 5, a = 4, — 0, and L = 4. As pp increases and fid decreases, 
the CDF of a Matern cluster process gets closer to that of a homogeneous Poisson point process with density 1 x lO"'^. 

for a, (7^, and Rd as above, a fixed number of antennas L = 4, and distance between the representative receiver 
and transmitter rx = 5, we consider Matern cluster processes with varying pairs of parent point density, pp, and 
mean number of daughter points /i^, per cluster. The pairs {pp,fid) are : (10^^, 1), (10~^, 100), (4 x 10~^,250), 
(2 X 10~^,500), and (10~^,1000) such that the effective density of interferers Pp pd in the network is fixed at 
10"^. Comparing the resulting CDFs which are shown in Figure |9j note that as the density of the parent points 
increases and the density of the daughter points decreases correspondingly, the CDF of the SIR in a Matern cluster 
network becomes similar to that of a homogeneous Poisson point process. This property can be explained by the 
fact that as the mean number of points per cluster in a Neyman-Scott process becomes small, the resulting process 
approaches a homogenous PPP because clusters with more than one daughter points become statistically negligible, 
and the interferers from clusters with single daughter points form a homogenous PPP by the displacement theorem 
for PPPs (e.g. see |161). 

VIII. Summary and Conclusions 

A technique to compute the CDF of the SINR on a link with a multiantenna receiver in a non-homogeneous 
Poisson field of interferers is presented and used to find expressions for the CDFs of the SINR for several 
representative non-homogeneous interferer distributions including power-law intensity functions for which a closed 
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form expression for the CDF is found, and circularly symmetric Gaussian intensity functions, for which the CDF 
is found for specific values of path loss exponents and is well-approximated by polynomial intensity functions. 
These results can be used to characterize the SINR in the center of a deterministic cluster. They also indicate that 
that while the SINR is significantly smaller in the center of a dense cluster than in a homogeneous network, the 
performance loss can be mitigated by using a larger number of antennas. 

We also used these results to show that if the number of receiver antennas is scaled linearly with the nominal 
interferer density in non-homogenous networks, the SIR converges in probability to a positive constant, indicating 
that it is possible to scale such networks by increasing the number of antennas, provided that the system assumptions 
hold. This particular finding generalizes a similar result for homogenous Poisson networks found in IS, Q and 
0. 

Furthermore, we explore Poisson cluster networks where clusters of interferers are randomly distributed, and 
provide theoretical CDFs for the SIR. These CDFs could be used to predict the performance of networks with 
clustered users where the cluster centers are not known a priori. 

The results in this paper generalize known results for homogenous Poisson networks with multi-antenna nodes to 
non-homogenous Poisson and Poisson cluster networks which are more realistic models for spatial node distributions 
in many cases. 



Appendix 

A. Proof of Theorem [7] 

Established by equations (11) and (12) in |6|, the CDF of 7 can be expressed as 



F^(7) = 1 - exp(-a27)E„ 
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where E^; represents expectation with respect to the random variable x. Recall that p = r " and the locations of 
the interferers are characterized by the PDF fr,e{r,6). Consequently, we have 
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As the number of interferers n is a mean Poisson random variable, the CDF of 7 is: 
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(43) 



Applying a sequence of steps similar to that used in the proof of the main result for homogenous networks in ll6l 
yields: 
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Given the relationship between the PDF of the locations of the interferers and the intensity function in (O, we 
denote -0(7) as: 

roo /•2tt 

(45) 
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Substituting ( |45l ) into (l44l) and applying the binomial theorem and the series expansion of the exponential function, 
yields: 



kUi-kV. 

= 1-2^ ^ exp(-V'(7) - o- 7) = 1 — — . (46) 

The last equality follows from equation (6.5.13) in |[20l . Finally, differentiating the CDF yields the PDF of 7. 
B. Proof of LemmaUl 

Consider following random variables U = X]fc=o ^ ~ Ylik=i ^k, where Vk are i.i.d. Poisson random 

variables with mean q. Note that 

Q(L, qL) = = Pr (t/ < L - 1) = Pr (f7 + Vo/{L - 1) < l) . (47) 

where the previous equality holds because Q{L,qL) is the probability that a Poisson random variable with mean 
qL is less than or equal to L — 1. By the weak law of large numbers as L — > 00, both C7 — )• g and Vo/(L — 1) — )• 
in probability implying that U + Vo/{L — 1) ^ q in probability. The latter implies that 

limPrff7+J^<nU/°' ^ ^ (48) 

L^oo \ L-1 ) 1, if M > q. 



Setting n = 1 in (1481 ) and substituting ( |47]) into the resulting equation completes the proof 

C. Proof of Theorem |2] 

Assuming that the noise is negligible, the CDF in ^ from Theorem [T] can be expressed as: 

r(L) 

Let q = iil^dj) in Lemma [T] Then, we have: 



j,^^^^^^_ty^^:^g^_ (49) 
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Additionally, since convergence in distribution to a constant implies convergence in probability (e.g. see Ell ). 7 
converges in probability as well. 

D. Proof of Lemma^ 

It is shown in Section 3.2 of |[T6l that given any positive measurable function ({.), 

En[eQ=] = exp I / p (e®e(-) _ l) du] , (51) 



for any negative and real ©. Taking the kth derivative with respect to G, (ISTb becomes 

^ mi!l!'^im2!2!'^2...mfc!A;!'"'= Vir^ / ' 

(mi,...,m,fc)eA^fc j=l 

The last step follows from Faa di Bruno's formula ll23l : 

5e^^(«(e))= E „.iii^.„.,i2!:»....„,!.i~ ^""""'""fa'Q'>-n(^'"(9) 

(mi,...,mfc)eXfc j=l 

where /(.) and g{.) are arbitrary measurable functions. Evaluating (l52l ) at = — 1 yields ([T3] ). 



E. Proof of Theorem \3\ 

According to Theorem [T] the CDF of the distance-normalized SINR 7, conditioned on IT is: 
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k=o ' \ x,en y 



(53) 



As we remove the condition with Fy{'y) = En[i^7(7| H)], the CDF of 7 becomes: 
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To evaluate this expectation with respect to the process the parent points 11, we assume zero noise power cj^ = 
and apply Lemma |2] with H = en C(-'^i)' ^^d thus ( [54l ) becomes 
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Converting the integrals from Cartesian coordinates to polar coordinates yields dTTl l. 
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General Power Law Term 
Lemma 3: Given the assumption that a > 2, 
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Proof: Using symbolic integration software, we can express the integral in terms of a Gauss's Hypergeometric 
Function: 
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Furthermore, we apply Euler's hypergeometric transformation EOl so that 
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Taking i?' — )■ 00 , we can simplify the equation given the last argument of the hypergeometric function is 1 |[20l . 
so that 
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The last step follows from the Reflection Formula of Gamma functions EOl . 



G. Proof of Theorem |?] 

Given any 5 > 0, let 5i be the ratio of 5 and the maximum value of r^^-p^. Then, since 
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According to the Stone-Weierstrass Theorem |[T9l . there exists A^{r,6) such that for all r in [0, i?], there is an 
integer N such that m > N implies 
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Then, for any 6 > 0, there is an integer N such that m > N implies 
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which shows that ^^^_2^ rA^(r, ^) is uniformly convergent to ^^^_2^ rA(r, 0) in {R,oo). Consequently, we can 
interchange the order of the integration and the limit, and therefore 
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Moreover, we express ([T9l ) as F^{'y) = /(V'm(7))» ^^d ^(7) = fiip{'y)). Since /(V'm(7)) ^ continuous function 
of V'm(7)> as m ^ 00, ^^(7) ^ V'(7)> implying 
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which completes the proof. 
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